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ABSTRACT 

A population of early star-forming galaxies is the leading candidate for the re-ionization of the uni¬ 
verse. It is still unclear, however, what conditions and physical processes would enable a significant 
fraction of the ionizing (Lyman continuum) photons to escape from these gas-rich galaxies. In this 
paper we present the results of the analysis of HST COS far-UV spectroscopy plus ancillary multi¬ 
waveband data of a sample of 22 low-redshift galaxies that are good analogs to typical star-forming 
galaxies at high-redshift. We measure three parameters that provide indirect evidence of the escape 
of ionizing radiation (leakiness): (1) the residual intensity in the cores of saturated interstellar low- 
ionization absorption lines, which indicates incomplete covering by that gas in the galaxy. (2) The 
relative amount of blue-shifted Lyman a line emission, which can indicate the existence of holes in 
the neutral hydrogen on the front-side of the galaxy outflow, and (3) the relative weakness of the [SII] 
optical emission lines that trace matter-bounded HII regions. We show that our residual intensity 
measures are only negligibly affected by infilling from resonance emission lines. We find all three 
diagnostics agree well with one-another. We use these diagnostics to rank-order our sample in terms 
of likely leakiness, noting that a direct measure of escaping Lyman continuum has recently been made 
for one of the leakiest members of our sample. We then examine the correlations between our ranking 
and other proposed diagnostics of leakiness. We find a good correlation with the equivalent width of 
the Lyman a emission line, but no significant correlations with either the flux ratio of the [OIII]/[OII] 
emission lines or the ratio of star-format ion rates derived from the (dust-corrected) far-UV and Ha lu¬ 
minosities. Turning to galaxy properties, we find the strongest correlations with leakiness are with 
the compactness of the star-forming region (SFR/area) and the speed of the galactic outflow. This 
suggests that extreme feedback- a high intensity of ionizing radiation and strong pressure from both 
radiation and a hot galactic wind- combines to create significant holes in the neutral gas. These results 
not only shed new light on the physical mechanisms that can allow ionizing radiation to escape from 
intensely star-forming galaxies, they also provide indirect observational indicators that can be used 
at high-redshift where direct measurements of escaping Lyman continuum radiation are impossible. 
Keywords: galaxies: evolution; galaxies: ISM; ISM: jets and outflows; 


1. INTRODUCTION 

Current evidence from observations of the cosmic mi- 


crowave background (CMB ; Dunkley et al. 2009 Planck 


Collaboration et al. 2014) a nd Lyman a fores t absorp- 


tion m quasar spectra (Ag. Fan et al. 2006) suggests 


that the reionization of the universe occurred over an 
extended period between redshifts of ^ 15 to 6. Deep 
near-IR (NIR) imaging with the Hubble Space Telescope 
(HST) indicates that the UV luminosity density of early 
star-forming galaxies is high enough that they are the 


photons necessary for this process (e.g. 

Bouwens et al. 

2012 

Shull et al. 2012 pesch et al. 2013 1 

kobertson et al. 

2015 
by w 

). For this to be possible, a mechanism is required 
nich Lyman continuum photons, produced by young 


stars in the centers of gas-rich star forming galaxies, can 
escape into the intergalactic medium (IGM) and the es¬ 
cape fraction of Lyman continuum photons fr om these 
galaxies must be relatively high: > 0.2 (e.g.|Bouwens| 

et al .112011] |Robertson et al. 2013). 
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Studies of high redshitt galaxies have yet to uncover 
a large sample of galaxies with the required high escape 
fractio n of ionizing photons. For example, Iwata et al.| 
(2009) only find evidence of escaping continuum photons 
in 17 of 198 high redshift {z ~ 3.1) Lyma n Break galaxies 
(LBGs) and Lyman a emitters (LAEs), |Vanzella et al.| 
(2010 2012) find evidence for escaping Lyman contin- 
uum photons in only on e galaxy of 102 LBG s studied at 
a redshift of z ^ 4 and Nestor et al. (2013) detect Ly¬ 
man continuum in 6 of 41 LBGs and 17 of 91 LAEs at 
z > 3.0. Even when escaping Lyman continuum photons 
are detected, in most high redshift studies /esc is usu- 
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et al. 12002 Fernandez-Soto et al.|2003 

Inoue et al.|2U05 

Shapley et al.||2006 Mostardi et ai.||201 

S and references 

above). This is higher than the upper limit (e.g. Lei- 

therer et al.|1995l IHeckman et al.|2001 

Deharveng et al. 

2001 Crimes et ai.|2009|) or value (e.g. 

Leitet et al.|2011 

20131 ot a tew percent measured in most local galaxies 

or our own Milky Way (Putman et al. 

2003) as well as 

the value measured in galaxies at redshitt z ^ 1 (e.g. 

Malkan et al. |2003| ISiana et al.||20071 |Cowie et al.||20091 


reionize the universe. Such a small number of detections 
is perhaps not surprising as an Hi column density of only 
1.4 X 10^^cm“^ at the Lyman edge is enough to produce 
unit optical depth and most star-forming regions from 
galactic disks to local starburst and high redshift star 
forming galaxies are gas-rich with column densities of 
10^^ — 10^"^ cm“^ ( Keimicutt||l998 Genzel et al.j|2010 ). 


One possible solution is to assume that galaxies have 
a “patchy” coverage of high column density, optically 
thick neutral gas clouds located in an optically-thin neu¬ 
tral medium. Such a “picket fence” model would allow 
ionizing photons to escape along optically thin lines of 
sight not blocked by an optically thick neutral cloud. 
In such a model, the escape fraction, /esc, of ionizing 
radiation is simply the percent of the galaxy not cov¬ 
ered by neutral clouds, fc (/esc = 1 ~ /c)- However, the 


relatively unusual (e.g. |Dove et al. 2000 

Clarke & Oey 

20021 IFuiita et al.U2003llWise & Cen||20Uy 

IHazoumov 

Sommer-Larsen||20 lOp . 


inherently difficult. At the epoch of reionization, the 
spectral region shortward of Lyman a is opaque making 
it impossible to study. Even for galaxies at a redshift 
above z ~ 4 the distribution of intergalactic neutral hy- 
dr ogen causes an atte nuation of > 1.8 mag below Lyman 


a (Inoue et al. 2014). In addition, studies of high red- 
shiit star-forming galaxies, with the exception of small 
samples of lensed galaxies, suffer from low spectral res¬ 
olution and poor signal/noise for individual objects. A 
sample of low redshift galaxies similar to high redshift 
star forming galaxies would allow us to search directly for 
escaping Lyman continuum photons with high sensitiv¬ 
ity, to calibrate/validate indirect indicators of escaping 
ionizing photons that could be employed at high red¬ 
shift, and provide clues as to the physical conditions and 
mechanisms that fo ster the escape of ioni zing radiation. 
A previous study (Heckman et al. ra hereafter Hll) 


used HST Cycle 17 data from the Cosmic Origins Spec¬ 
trograph (COS) in the far-UV (FUV; HST Program 
11727: PI T. Heckman) to study 8 starbursting local 
galaxies known as Lyman Break Analaogs (LBAs). Due 
to the relatively poor sensitivity of COS beyond the Ly¬ 
man edge for objects with z ^ 0.1 — 0.25, these obser¬ 
vations were made at longer wavelengths and the escape 
fraction of ionizing radiation wa s inferred from s everal 
indirect indicators. Previously, Overzier et al. (2009 


hereafter 009) and Hll compiled a list of indirect indi- 
cators that may provide evidence for the presence of es¬ 
caping ionizing continuum photons. These include: 1) a 
large ratio of the star-formation rate (SFR) derived from 
dust-corrected UV to that derived from dust-corrected 
Ha possibly due to the loss of ionizing photons. 2) weak 


[S II] emission relative to Ha or [N II], possibly indicat¬ 
ing matter-bounded HII regions. 3) low-ionization in¬ 
terstellar absorption lines that were optically-thick but 
with significant residual intensity in their cores, indica¬ 
tive of “holes” in that gas. And 4) Lyman a emission 
lines with a significant ratio of blue-shifted:red-shifted 
emission again indicative of holes in the neutral ISM that 
would permit leakage of both Lyman a emission and ion¬ 
izing continuum radiation. 

Hll also found that objects with a dominant central 
object (DCO; 009)- a massive (M* > 10®Mq) and 
very compact (R^lOOpc) starburst- were more likely to 
show strong supernovae-driven winds (~1000 kms/s) and 
“patchy” neutral hydrogen coverage with a range of dust- 
free escape fractions from 0.3 < ^ 0-6- How¬ 

ever, this sample was not large enough for a full statis¬ 
tical analysis. Thus, in HST Cycle 20 we successfully 
proposed for COS FUV-spectroscopy of an additional 14 
LBAs (HST Program 13017 : PI T. Heckman). In this 
paper we report the results of an analysis of these new 
data and re-analysis of the previous Cycle 17 data which 
allows us to draw statistically robust conclusions about 
the degree to which these different indirect indicators of 
escaping ionizing radiation agree with one another and 
ascertain which properties of the starburst correlate best 
with these indicators. 

Our results are made even mor e relev ant by re¬ 
cent results from fBorthakur et al. (2014) that show 
direct evidence for escaping Lyman continuum from 
one of the LBAs in Hll. The object in question, 
SDSSJ092159.38-k450912.3 (J0921) contains a DCO with 
a relatively high (~ 25%) predicted escape fraction of 
Lyman continuum photons from the variety of indirect 
methods listed above. Follow-up observations using COS 
below the Lyman edge found a luminosity, XL\ at a rest 
wavelength of 910A of 5.0 x 10"^^ erg s“^, implying a 
dust-free escape fraction of = 21 ± 5%, similar 

to that predicted by indirect indicators. 

The rest of the paper proceeds as follows; Section § 
describes our sample selection, § and § describe 
our observations and data analysis. Section § com¬ 
pares the indirect indicators of escaping Lyman contin¬ 
uum photons with each other and §|^ examines their de¬ 
pendence on the galaxy/starburst properties. We offer 
our conclusions in Section § The names of all ob¬ 
jects in the paper are derived from the hour and min¬ 
utes of their J2000 right ascension coordinates as taken 
from the SDSS. Full coordinates are provided in Table 
1. Throughout we adopt a flat ACDM cosmology us- 
ing Da = 0-7, Dm = 0.3,and Hq — 70 km s“^ Mpc“^ 
( Spergel et al.|[2007 ). 


2. SAMPLE SELECTION 


As described above, the existence of a local popula¬ 
tion of LBA galaxies has be en established over the past 
decade (Heckman et al.|2005 ). This sample was originally 
selected by combining optical data from the Sloan Dig¬ 
ital Sky Survey (SDSS; [York et ^ 2000) with UV data 


from the G alaxy Evolution Explorer (CALEX; Martin 
et al. 20051 to match the UV properties of LBGs using 
the sample criteria: 

1. No detectable AGN contribution to continuum 
emission 























































































Indirect Evidence for Escaping Ionizing Photons in Local Lyman Break Galaxy Analogs 


3 


2. FUV luminosity of Lpuv > 2 x I0 ^°Lq measured 
at 1550A 

3. FUV effective surface brightness (mean FUV sur¬ 
face brightness within the half-light radius of the 
SDSS u-band) of Ipuv > 10®Lq kpc ^ 

Subsequent research has shown that this LBA sam¬ 
ple is matched with LBGs in terms of morphology, size, 
UV luminosity, SFR, mass, gas frac tion, velocity disper¬ 
sion, metallicity, and dust content (Hoopes et al. 2007[ 
'Basu-Zych et al.||2007[ [Overzier et al.||2UU8t |Basu-Zych| 


et al. 1120091 lOverzier 

et al.||2009| 

Overzier et al.||2010 

2011 Gonga 


Goncalves et al.|^10[ 

longaives et al.|2014p. Thus, 

ok tor clues as 
to how outflows of gas and ionizing radiation from early 
star-forming galaxies enriched and reionized the IGM. A 
list of all the LBAs observed in our two HST campaigns 
with COS along with details of the observations and their 
properties is provided in tablef^ Objects were chosen for 
follow-up with COS that had the highest FUV continuum 
fluxes (as evidenced by their GALEX photometry) and 
appeared to be the most compact (this applies only for 
the sample from Hll where previous UV imaging from 
HST was available) so as to provide the best possible 
spectral resolution. The sample presented here is based 
on LB As selected in the original sample of [Heckman et ah] 


(2005), supplemented with new targets found in the lat- 
est crossmatch (Ove rzier et al. in prep.) b etween SDSS 


data release 7 (DR7;|Abazajian et al.|2009) and GALEX 
GR6. 


3. HST OBSERVATIONS AND DATA ANALYSIS 

Our observations were conducted using the COS in¬ 
strument aboard HST in cycle 20. Each observation in¬ 
cludes a 100 sec acquisition image in the near-UV(NUV) 
as well as a longer spectroscopic exposure (1-2 orbits 
depending on the FUV flux expected to pass through 
the COS aperture). Our spectral observations use the 
COS FUV G130M and G160M medium resolution grat¬ 
ings (R~17000). This spectral resolution ensures that 
the spectral line-spread function is smaller than the ve¬ 
locity dispersion of the ISM within the LBA and so al¬ 
lows us to resolve important absorption- and emission 
line features. Our observations span the observed wave¬ 
length range from A = 1160A — 1780A and are de¬ 
signed to allow maximum spectral coverage from rest- 
frame O VI (A 1035A) which probes coronal gas at 
T~ 10^ - 10® K to C IV (A 1550A) and so probes 
the O star population through stellar winds. In ad¬ 
dition, this wavelength range covers several important 
ISM transitions such as Lyman a, the low-ionization 
transitions of Si II (A 1190,1193,1260,1304,1526A) and 
C II (A I 334 A) and the intermediate ionization transition 
of Si III (A 1206A). Of the newest sample six objects 
are at z < 0.11 meaning the O VI line is out of reach 
while three are at z > 0.16 meaning C IV is not covered. 
Postage stamp images of the new sample of 14 LBAs in 
our HST COS Cycle 20 sample are shown in Figures 
and Note that the images shown are the acquisition 
images and that those in figure were taken with the 
MIRROR B setting and so show a secondary image on 
the detector. UV images of the previous Hll sample can 
be found in 009. 


3.1. COS Data 

All spectra were run through a suite of procedures sim¬ 
ilar to those described in Hll. Data were first retrieved 
from the MAST archive and were processed through the 
standard COS pipeline. Then, the G130M and G160M 
spectra were merged and rebinned by 21 pixels (corre¬ 
sponding to a spectral bin size of ^ 50km s“^); no at¬ 
tempt was made to smooth the data at this point. The 
data were then normalized by fitting a power law fol¬ 
lowed by a low-ord er polynomial function to t he contin¬ 
uum as defin ed by Calzetti et al. (19941 and Leitherer 
et al.j (|2002| defined m regions of the spectra designed 
to avoid strong emission and absorption features both in 
the observed and rest frame) and then dividing the spec¬ 
tra by this polynomial. The full sample was combined 
to make a single high signal-to-noise spectrum shown in 
figure 1^ 

Our goal in this paper is to measure the properties of 
the interstellar gas. To that end, we undertook a proce¬ 
dure to remove the spectral features produced by stars 
(e.g. photospheric a nd stellar wind lines). Starburst99 
(from now on SB99; Leitherer et al. 1999) was used to 
generate synthetic spectra based on stellar evolutionary 
synthesis models. We produced models based on two 
star formation histories: an instantaneous starburst with 
no subsequent star formation and a continuous starburst 
with a constant rate of star formation. The stellar popu¬ 
lation was parameterized by a Kroupa initial mass func¬ 
tion (IMF; iKroupa 2001). The metallicity of the stellar 
population was previously determined by 009 and we 
created two sets of model spectra with SB99 to bracket 
the observed metallicities. The stellar population then 
evolves from the zero-age main sequence using the high 
mass loss evolutionary models of the Geneva Group. The 
model spectrum we used are generated from model at¬ 
mospheres, not an e mpirical library; these spectra are 
described in detail in Leitherer et al. (2010|. We sam¬ 
ple the model every two million years to generate the 
synthetic spectra for our model fitting. Our observed 
COS data were interpolated to the wavelength grid of 
the spectra produced by our stellar evolutionary synthe¬ 
sis models. This caused some degradation in our spectral 
resolution, especially at shorter wavelengths. 

The method we e mployed for findin g the b est fit model 
is also described in Tremonti et al. (2001). First, we 
do our own normalization of the SB99 models accord¬ 
ing to the same prescription used for the COS data. 
We compare our data and models based only on a 
select few stellar wind lines generated by massive O 
stars. For this analysis we used the regions surround¬ 
ing C IIIA II 75 A (1171 - II 8 OA), N VA 1239A(1225 
- I 247 A), and C IVA 1548A (1534 -1560A), except for 
those objects for which C IV is not available. In two 
cases (J0823, J1113) a damped Lyman a (DLA) system 
and circum-galactic medium of a nearby galaxy respec¬ 
tively obscured the N V feature, preventing us from using 
it in our comparison. In one of these cases, J1113, there 
was also no C IV available as the object was at a higher 
redshift and so the model comparison was tied to a single 
absorption feature. There is a known problem with the 
O VI A 1038A feature in SB99 due to the Doppler broad¬ 
ening of the Ly/3 emission line generated by SB99 that 
prevents us from using it to fit our data. However, we do 
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Figure 1. NUV images of the new LB A sample. These images were taken using COS in 100 sec exposures primarily for 
acquisition purposes. Each image is 2” per side with the corresponding physical scale for 1” marked. The objects were imaged 
with COS in the MIRRORA configuration. 


find that our data shows rough consistency with the SB99 
O VI when this line is excluded from the fit. We then 
calculated the chi-squared value over all selected regions 
for each model spectra and take the “best-fit” model to 
be the model with the lowest chi-squared value. We find 
that, in most cases, our LBAs are best fit by models 
with ages less than 10 million years whether instanta¬ 
neous or continuous starbursts (see figure]^. We then 
subtract the best fit model from our data to remove the 
stellar spectral component. With the stellar continuum 
and line features removed, we are able to study the ISM 
signatures in detail. 

Unless otherwise stated all measurements can be as¬ 
sumed to have errors on the order of 10%-15% dominated 
by systematics in the polynomial fit to the continuum 


emission and subtraction of the SB99 models. 

3.2. Measurements of velocity outflows 

Strong outflows of gas are generic in intensely star¬ 
forming galaxies at low- and high-redshift. One signature 
of these outflows are the broad blueshifted interstellar 
absorption lines (see figure [^. Hll suggested that very 
high outflow speeds might oe associated with galaxies 
with escaping Lyman continuum emission. In order to 
characterize the outflow speeds we use the Silll A I 206 A 
transition as opposed to Si II or C II because it is typi¬ 
cally the strongest interstellar absorption feature in the 
spectra of the LBAs allowing us to characterize outflow 
speeds even in cases where significant residual intensity 
makes other lines shallow. In the case of JI414 where 
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Figure 2. NUV images of the new LB A sample. These images were taken using COS in 100 sec exposures primarily for 
acquisition purposes. Each image is 2” per side with the corresponding physical scale for 1” marked. The objects were imaged 
with COS in the MIRRORB configuration. MIRROR B produces a secondary image seen here to the right of the primary image 
for each object. 

Si III was not available in the spectrum, the outflow ve¬ 
locity was calculated using the C II absorption feature. 

We measure the flux-weighted line centroid of this fea¬ 
ture, integrating over a spectral region selected by eye 
(shown in Figure in grey) to have detectable absorp¬ 
tion. Future work (Heckman et al. 2015, in prep) will 
involve a more detailed analysis of the outflow velocities 
derived across the range of available ISM transitions. We 
found a maximum outflow speed of 561 km s“^ (J0808) 
with an average outflow speed of 375 km s“^. Infilling of 
absorption lines from scattered emission could cause an 
offset between the measured and actual outflow veloci¬ 


ties though t his effect is likely small (see section 3.4 of 
this paper or |Scarlata & Panagia||2015| |Prochaska et al.| 
^OTTl. 


3.3. Residual Flux in saturated ISM lines 

Hll showed that the strong absorption lines tracing 
the low-ionization ISM in LBAs are saturated (optically 
thick) but do not always have zero residual intensity in 
the line center. They interpreted this as evidence for par- 


neutral gas (as in e.g. Heckman et al. 

2000 Rupke et al. 

2002 Quider et al.|2009 

2010). In this model, optically- 


thin medium and thus the percentage of the ionizing 
source not covered by neutral gas clouds is the same as 
the residual flux (percent) in saturated ISM lines. This 
technique is increasingly used at higher redshift where 


permitted (see, for example Jones et al. 2013). 

This is most clearly seen by examining the series 
of Si II transitions which span over an order of mag¬ 
nitude in oscillator strength (f): A 1190.416A(f = 
0.277), A 1193.2901(0.575), A 1260.4221(1.22), 
A 1304.3701(0.093) A 1526.707l(0.133)b If the 
UV sources were being viewed through a uniform 
covering of optically-thin clouds we would expect the 
residual intensity in each line to decrease as the oscil¬ 
lator strength (optical depth) increases. On the other 
hand, if our UV sources were only partially covered by 
optically-thick clouds, there would be no dependence 
of residual intensity on oscillator strength and instead 
the residual intensity at the core of each line (/q) would 
be related to the covering factor (A) of neutral clouds 

{lo/Icont = 1 - /c). 

We began by fitting the 1190l, 11931, and 1260l 
lines, which should be the most optically thick. These 
absorption lines were fit with a Gaussian m odel using 
the IDL package MPFIT (Markwardt 2009). We con¬ 
strained all three lines in a tit to have tne same centroid 
and width in velocity space but allowed the line strength 
to vary so as to be able to study changes in the resid¬ 
ual flux across different absorption lines. Residual flux 
was measured for each line by calculating the depth of 
the line at the centroid calculated from our model. We 

^ oscillator strengths from NIST Atomic Database: 
http: //www.nist.gov/pml/data/asd.cfm 
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Figure 3. Coadd of our full sample of 22 COS spectra. No weighting of individual spectra was applied. All spectra have been moved 
to the same wavelength grid as the SB99 models. Important ISM lines are noted above in black and important stellar features are noted 
below in red. Lines with both stellar and interstellar components are indicated in both colors. 


found that the residual flux in all three lines is largely 
consistent, and frequently large. 

To further confirm that the strong Si II transitions 
in our spectrum were optically thick, we calculated the 
evolution in equivalent width (EW) for the five different 
Si II transitions (where available within the wavelength 
range of each individual spectrum) as a function of os¬ 
cillator strength (commonly referred to as the “curve of 
growth”). This allows us to search for the expected evo¬ 
lution from the optically thin to optically thick regime. 
We hnd for many galaxies where all five transitions are 
visible in the spectrum, there is a clear transition from 
the optically thin portion of the curve of growth where 
the equivalent width evolves linearly with A/, to the op¬ 
tically thick regime where the curve of growth flattens. 
This is shown for our high signal-to-noise LBA spectrum 
in figure Note that there is a transition of S III at 
1190.206A that contaminates the Si 11 1190.416A ab¬ 
sorption line and probably explains why it unexpectedly 
appears stronger than Si II 1193 in our figure though the 
latter has a stronger oscillator strength. This gives us 
confidence that for stronger absorption features of Si 11, 
(in particular, A = 1260) we can assume that the line 
is optically thick. Thus, the relative residual intensity 


in the core of the absorption line can be interpreted as 
a measure of the fraction of the UV source not covered 
by optically thick low-ionization gas. It should be noted 
of course that while in our model these sightlines are 
optically-thin to low-ionization metal lines, they could 
still be optically thick to ionizing radiation. Thus, ob¬ 
servations below the Lyman edge are required to confirm 
the actual fraction of escaping ionizing radiation. 


3.4. Resonance and Fine Structure Line Emission 


As emphasized by 
& Panagia (2015), t 


Prochaska et al. (20111 and Scarlata 
re U V absorption lines we are using 
are resonance lines and as such, each upward transition 
resulting in absorption should be followed by a downward 
radiative transition resulting in emission. This emission 
will ‘infill’ the absorption line, and can in principle have 
a significant effect on the net observed profile. Thus, cor¬ 
rectly measuring properties such as optical depth, cov¬ 
ering factor, and outflow velocity requires correcting the 
profiles for this infilling. 

Both of the above papers show that the fine-structure 
transitions that are associated with many of the 
commonly-used resonance lines can be used to deduce the 
properties of the infilling resonance line emission. These 
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J015028 Age = 4.0Myr Model = Kroupo IMF, Instantaneous Burst, Solor Metollicity 
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Figure 4. Best fit model for a single LBA, J0150. The data is shown in black and over plotted in colonr is the best fit model 
originally generated by SB99. The model has three components; in red is the model continuum (not used in the fit), in blue 
are regions of the model that are masked because they overlap with interstellar featnres (not used in the fit). Finally, in green 
are the stellar wind li nes we used for our chi-squared test. This figured was generated using code from Dr. Christy Tremonti 
( Tremonti et al.||2004 |. 


fine structure lines share the same upper level as the as¬ 
sociated resonance lines, but the lower level is slightly 
higher in energy than the ground state associated with 
the resonance line. Thus the downward radiative transi¬ 
tion can occur via either the resonance line (resulting 
in infilling) or via the fine-structure line (which does 
not contaminate the observed absorption line). Quan¬ 
titatively, the flux ratio between the resonance and fine- 
structure emission should be the ratio of the respective 
Einstein A coefficients. 

In Figure we show two pairs of resonance and fine- 
structure lines associated with Si II: Si II 1260.4 plus Si 
II* 1264.7 and Si II 1526.7 plus Si 11*1533.5 in the coad¬ 
ded spectra of our 22 LB As. The ratios of the Einstein 
As (resonance/fine-structure) for these two pairs are 0.87 
and 0.051 respectively. We measure respective ratios of 
the equivalent widths (resonance/fine-structure) of 0.10 
and 0.11, implying that the respective ratios of the in¬ 
ferred (infilling) emission to observed absorption is only 


0.09 and 0.006. The fine-structure emission lines are nar¬ 
row and centered on Vgys- It is therefore clear that in¬ 
filling of the absorption lines will have a negligible effect 
on the properties of the gas inferred from the observed 
absorption lines. 


3.5. Lyman a Emission 

Hll suggested that Lyman a emission line profiles 
in which a significant fraction of the emission was 
blueshifted relative to the galaxy systemic velocity were 
indicative of galaxies in which there were holes in the dis¬ 
tribution of outflowing neutral hydrogen. Indeed in our 
galaxies we notice a variety of different Lyman a profile 
shapes. Many galaxies display the traditional P-Cygni 
profile with redshifted emission and blue shifted absorp¬ 
tion t hat are commonly seen in LBGs (|Shapley et al. 
2003 Steidel et al. 2010). However, other LBAs dis¬ 
play prohles with an abundance of blueshifted Lyman 
a emission. Following Hll, we have quantified the dif- 
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Figure 5. Spectral region containing Si IIA 1190,1193A and Si IIIA 1206A absorption features tracing the low- and mid¬ 
ionization ISM. The spectra shown have been normalized and the best-fit stellar model subtracted so the normalized continuum 
flux is equal to 1. Areas where MW features interfere with the spectra have been masked. The rest wavelength for each line, 
determined by the SDSS redshift, is marked in black. The gaussian model fit to both Si II transitions is shown in red. Each 
absorption line is fit with a single gaussian component constrained such that both transitions as well as Si IIA 1260A have the 
same width and centroid but the flux in each transition is allowed to differ. The region over which the outflow speed is calculated 
from Si III is indicated in grey. In the case of J1414, where Silll was not available in the spectrum, the outflow velocity was 
calculated using the CII absorption feature. Notice the significant amount of residual flux in both Si II lines as well as the 
asymmetric, blue profiles of all three lines. 
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Figure 6. Left: Evolution of Si II transition rest-frame equiv¬ 
alent width as a function of oscillator strength (log A/) also 
referred to as the “curve of growth” for the co-added spectrum 
of our 22 LBAs. We see the characteristic evolution from opti¬ 
cally thin absorption (EW/\ oc /A) to optically thick absorp¬ 
tion where the curve-of growth becomes fiat. This provides 
confidence that our strongest Si II absorption lines, particu¬ 
larly A 1260 is indeed optically thick and any residual inten¬ 
sity in the core of the feature must be due to a less than 100% 
covering of neutral gas. Note that there is a transition of S III 
at 1190.206A that contaminates the Si II 1190.416A absorp¬ 
tion line and probably explains why it unexpectedly appears 
strong than Si III 193 in our figure. Right: Model fit to the 
three most prominent Si II transitions. Our model uses a sin¬ 
gle gaussian for each absorption line constrained such that all 
lines have the same width and centroid, only the flux in each 
transition is allowed to differ. 




Rest Wovelength (Angstroms) 


Figure 7. Zoom-in on the spectral regions around Si II 1260.4 
& Si IP 1264.7 (left) and Si II 1526.7 & Si IP1533.5 (right) 
from a combined spectra of our 22 LBAs. Both absorption 
lines have associated resonance and fine-structure lines that 
can be used to estimate the amount of in-filling expected 
in the absorption line. We measure respective ratios of the 
equivalent widths (resonance/fine-structure) of 0.10 and 0.11, 
implying that the respective ratios of the inferred (infilling) 
emission to observed absorption is only 0.09 and 0.006 and 
thus has a negligible effect on the properties of the gas in¬ 
ferred from the observed absorption lines. We show in grey 
on the right where a local continuum was used to measure 
the equivalent width of Si 11*1533.5 otherwise the inferred 
value of infilling emission would have been underestimated. 


ferent Lyman a profiles shapes in the following way: we 
first measure the Lyman a equivalent width blue-ward 
and red-ward of the systematic velocity, where we de¬ 
fine the equivalent width for emission (absorption) to 
be positive (negative). We then take the ratio of the 


blue/red equivalent widths to parameterize the profile 
shape (Reqw)- Therefore, a negative value of Reqw indi¬ 
cates blue shifted absorption and redshifted emission (a 
traditional P-Cygni profile) while large, positive values of 
Reqw indicate strong redshifted and blue shifted emission 
(see the comparison in figure [^. It is worth noting that 
all but one of our objects shows strong Lyman a emis¬ 
sion, a s opposed to some hi gh redshift samples such as 
that of|Shapley et al.| (|2003|). 

The Lyman a equivalent width itself is of interest since 
it is a measure of the fraction of Lyman a photons that 
escape the galaxy relative to the escaping fraction of 
FUV continuum photons. It is therefore plausible that 
the equivalent width of the Lyman a emission line could 
be used as an indirect indicator of the escape of Lyman 
continuum photons. We have therefore measured a net 
Lyman a equivalent width in emission for comparison 
with high-redshift samples. 

For several objects (J0055, J1113, J1144, J1414, and 
J1416) the continuum normalization around Lyman 
a was poor due to intervening absorption from either 
the circumgalactic medium of nearby galaxies or the lo¬ 
cal interstellar medium. J1113 clearly shows evidence 
for a DLA at z = 0.1712 associated with the CGM of 
J111324-1-293051 with metal lines that are all double- 
peaked and a strongly blue-shifted component that is 
due to the foreground galaxy also seen in the SDSS im¬ 
age at the same redshift. In the other four cases, that of 
J0055, J1144, J1414, and J1416, we see only very weak 
Lyman a emission-lines. One possibility is that this is a 
result of errors in the global continuum fitting because 
there is no true continuum around Lyman a. Figure 1^ 
demonstrates that the spectrum around Lyman a ±25A 
or so could be affected by the broad N V stellar wind line 
redward and the blending of absorption due to Lyman a, 
Si III and N I on the blue side. Alternatively, this pro¬ 
file could also be due to underlying Lyman a absorption 
in the same galax y as seen in, for example, [Schaerer fc| 
Verhamme (20081 and Dessauges-Zavadsky et al. (|2010p. 


In either case, we choose to calculate Reqw and the 
equivalent width with a local continuum estimated in a 
region close to the Lyman a emission line and fluxes were 
calculated relative to this local value. If the broad Lyman 
a absorption is produced by a foreground galaxy (exter¬ 
nal to the outflow) or by the dense interstellar medium 
of the galaxy (internal to the outflow), then we believe 
that measurements of the outflow component of Lyman 
a should be made relative to the local continuum allow¬ 
ing the measurement to highlight the portion of Lyman 
OL emission that is escaping from the galaxy. In Figures 
[Ml] and we show the objects where a local con¬ 
tinuum was used in purple to highlight that they are 
not driving the observed correlations. Measurements of 
Reqw using a global continuum could not be made in a 
meaningful way for comparison as it is not clear over 
what range the intervening absorption makes a signifi¬ 
cant contribution to the spectrum. We find that these 
five measurements do not significantly affect our conclu¬ 
sions. 


3.6. Starburst Size and Morphology 

All morphological information for our new sample of 
LBAs is taken from acquisition images for COS which 
are performed using the NUV MAMA detector. The UV 
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Figure 8. Detailed Lyman a profile of all LBAs in this sample. Zero velocity is marked with a dashed red line. The area over 
which values of Reqw and EW were determined are marked in grey. These were chosen by eye. For well-normalized data the 
continuum is expected to have a value around 1.0 (grey dotted line). In the cases of J0055, J1113, J1144, J1414, and J1416 the 
global continuum fitting produced a poor fit around Lyman a so a local continuum was estimated in a region close too Lyman 
a and Req-w was calculated relative to this local value. There is significant variety in the observed profile shapes. While many 
objects display the traditional P-Cygni profile expected (resulting in a negative or small positive value of Reqw), several show 
both redshifted and blueshifted emission resulting in a large positive value of Reqw- 
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half-light radius was determined using SExtractor to es¬ 
timate the half-light radius from the Kron radius. The 
measured half-light radii are listed in table For the 
previous sample UV half-light radii are measured from 
existing HST data taken with the Solar Blind Channel 
(SBC) on the Advanced Camera for Surveys (ACS) of 
HST (see 009 section 2 for more details) using SExtrac¬ 
tor. 

Most LBAs consist of a number of compact star¬ 
forming clumps with a range of sizes and brightnesses. 
However, a minority consist of only a single, highly com¬ 
pact (^100 pc) and massive (~ 1 x 10®Mq) central 
source. Hll and 009 have dehned these “dominant cen¬ 
tral objects” (DCOs) as those LBAs with extremely com¬ 
pact star bursting central regions as determined by UV 
photometry. Any object where the ratio of UV flux in 
the central 0.2” to the flux within the FWHM of the 
COS PSF is greater than ^ 25% is considered a DCO 
(009; see Figureand [^. This calculation was done for 
our Cycle 20 data using the Aperture Photometry Tool 
(Laher et al.|2012|) on the HST acquisition images to cal- 
cuiate the flux within each radius. The results of this 
calculation can be found in table [21 The classifications 
for the previous sample were taken from 009 and con¬ 
firmed both visually and using APT. Hll found a connec¬ 
tion between those LBAs that were also DCOs and those 
objects that were strong candidates for escaping Lyman 
continuum photons hypothesizing that strong starbursts 
were capable of creating optically thin channels in the 
neutral ISM that would allow ionizing photons to escape 
the galaxy. 


4. ancillary data 

4.1. Star Formation Rates 


It was originally hypothesized in 009 that a system¬ 
atically smaller SFR as calculated from Ha compared to 
SFR indicators from the FUV and infrared (IR) might 
indicate that ionizing radiation was escaping the galaxy. 
While the FUV and IR luminosities together should be 
a robust tracer of the bolometric luminosity and SFR, 
escaping ionizing radiation might depress the measured 
value of the extinction-corrected Ha luminosity. 

FUV luminosities are measured at rest-frame I500A 
from our COS data and then converted to 1600A (not 


available in our dataset) assuming f\ oc where /3 is the 
spectral slope {(3 = 2.321ogio(mi?c/y — mjsfuv)-'^-^) and 
the FUV and NUV magnitudes are taken from GALEX. 
This is a reasonable approximation in the FUV where the 
slope is expected to be roughly constant. The FUV lu¬ 


minosity is then simply Lpjjv = ALa at 1600A (Overzier 

eral][2^ . 


h'or the Hll s ample the IR luminosities are taken from 


Overzier et al. (2010) which uses SED fitting of infrared 
data obtained with the Infrared Array Camera (IRAC), 
the Multi-band Imaging Photometer (MIPS) and the In- 
fraRed Spectrograph (IRS) on the Spitzer Space Tele¬ 
scope to estimate the IR luminosity, Liji{3 — lOOO/rm). 
Unfortunately, none of the new sample in the current 
paper have been observed with Spitzer or Herschel, so 
our only infrared data comes from the Wide-Field In¬ 
frared Survey Explorer (WISE) All-Sky Survey where the 
W4 band is at an observed wavelength of 22fj,m. This 
data is dominated by hot dust in the mid-IR whereas 


the IR luminosity is typically dominated by colder dust 
seen at longer wavelengths. We have therefore deter¬ 
mined an empirical relation between IR luminosity mea¬ 
surements based on Spitzer plus Herschel and measured 
W4 luminosities for the 28 available LBAs (Overzier et 
al. in preparation). We use this relation to estimate 
the total IR luminosities in our new sample. Finally, the 
S FRpjjVrnrr is Calculated a ccording to the prescription 
of Kennicutt & Evans (2012) by first correcting the FUV 
luminosity with the IR lurninosity and then converting 
to a SFR. 


log SFRpjjv, corr {Mq yr 1) = 

log {Lppv,obs + 0.46 X L7fl,3_iooo/im( i^Aj,)) — 43.35 

( 1 ) 


To measure the Ha luminosity the emission must Hrst 
be corrected for dust exti nction. This is accom plished us¬ 
ing the reddening law of Calzetti et al. (2000) for wave¬ 
lengths greater than 6300A to generate the correction 
factor,‘f’ : 


/ = 10 o(°-4^^(®“^V^-®^®(“1-857+10400/A)+4.05) ^ 2) 

where the extinction, E{B — V) is calculated from 
the ratio of Ha to H/3 flux as E{B — V)gas = 

l°glo[(/Ha//ff/3)/2.87] 

0.4x1.17 

To convert our extinction-corrected Ha luminosity to 
a SFR we use the same prescription as 009: 

SFRho^AMo y^~^\ = 5-3 X s”'] (3) 

As in 009, an additional factor of ~ 1.7 is applied 
to account for the Ha flux outside the SDSS fiber. We 
use the Hnal indicator, log as possible ev¬ 

idence of escaping ionizing radiation. 

4.2. Measuring Stellar Masses 

We use both a stellar mass for the entire galaxy and 
for the starburst. The total stellar mass is taken from 
the MPA-JHU galaxy catalog of SDSS DR7^. We fit 
a burst mass by normalizing the starburst mass of an 
instantaneous burst (or SFR yr“^x age for a continuous 
burst) of our best-fit SB99 model for each galaxy to the 
flux from the UV spectrum. 

4.3. Measuring Emission Line Diagnostics 

It was previously noted in 009 that some LBAs, 
were characterized by unusually weak [S II]A6717,6731 
line emission compared to typical star-forming galaxies. 
Since in HII regions this line arises in gas at the outer 
boundary of the Stromgren sphere where the gas transi¬ 
tions from ionized to neutral (e.g. Pellegrini et al. 2012), 
relatively weak [SH] emission might indicate the presence 
of gas that is optically thin to ionizing radiation (e.g. 
density-bounded conditions). 

To quantify the relative strength of [S II] emission 
in the LBAs, we first defined the locus of normal 
star-forming galaxies in the diagnostic diagram of log 

^ Available from: http://www.mpa- 

garching.mpg.de/SDSS/DR7/ 
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Figure 9. Emission line ratio of log ([S II]/Hq:) vs. log([OlIl1 /17/3). The density c ontours are SDSS galaxies categorized as star 
forming galaxies according to the emission line diag nost ics oflKe-wley et al. (20061 in the BPT diagram. The dashed line is our 
empirically determined “star-forming” ridgeline (see |4.3[ |. We calculated the perpendicular distance between our objects, shown 
in red, and the mode of star-forming galaxies on the diagram to express the deviation from a typical star-forming galaxy for this 
emission line diagnostic. J0921, our object with confirmed Lyman continuum emission, is represented as a yellow diamond. We 
also include some other s amples for comparison. The other LBAs observed with HST in 009 are shown in black. We include 


the hig h-redshift SMG of|N esvadba et al. (|2007| in purple as well as the high-redshift star-forming galax ies of the MOSDEF 
survey (Shapley et al.||2014|) in green. F inally, "we include a sample of nearby warm IR-luminous galaxies I Kewley et al.||2001 ) 
in gray. 


([S II]/Ha) vs. log([Oiii]/i7/3). We took only those 
galaxies from SDSS that fell in the star-forming r egion of 


the log ([N II]/Ha) BPT diagram as defined by Kewley 
et al. (20061). We then followed a similar prescription to 


Kewley et ah (2006) to define the ridge-line. We began 
with an empirical end-point at the bottom of the distri¬ 
bution of SDSS “star-forming” galaxies and divided the 
sample in to radial bins of 0.1 dex. For each radial bin we 
calculated the mode of the angle between the end-point 
and each galaxy in the bin. This provided a parame¬ 
terization of the ridge line as shown in Figure The 
displacement value. Dp is then defined to be the perpen¬ 
dicular distance between each LBA and the parametric 
“ridge line”. It is evident in this figure that a number of 
LBAs tend to be significantly displaced from the ridge of 
normal SDSS star-forming galaxies. 


5. RESULTS 

5.1. Comparing Indirect Indicators of Lyman 
Continuum Escape 

Following 009 and Hll, we consider three parameters 
as being indirect indicators of escaping Lyman contin¬ 
uum photons: 

1. The amount of residual intensity in the core of the 
Si II A1260 ISM absorption line (%) 

2. The shape of the Lyman a line emission line profile 
(^Reqw ) • 

3. The distance from the star-forming ridge line in the 
space of log ([S II]/Ha) vs. log([Oiii]/iL/3) (dex) 

The calculated values for each indicator measured in 
every LBA from our sample are presented in table 
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When we compare each of these indicators, we find sta- 
tistically significant correlations between all of them (see 
Figure 101. We see a regime with high residual flux (indi¬ 


cating incomplete coverage of the starburst by optically- 
thick low-ionization gas), the presence of a significant 
amount of blue-shifted Lyman a photons escaping the 
galaxy (indicating incomplete coverage of the near-side 
of the outflow by neutral hydrogen), and a relative deficit 
of [S II] emission (indicating the presence of portions of 
the ionized gas that are optically-thin to the Lyman con¬ 
tinuum). As previously noted by Hll, these indicators 
tend to be most pronounced in those objects classified 
as DCOs with strong compact star forming cores. We 
will reassess this quantitatively in the next section. We 
do not find a strong correlation between these indicators 
and the ratio of SFRs derived from the FUV luminos¬ 


ity vs. the Ha luminosity (see section 5.21 and so have 
excluded it here as an indicator of escaping Lyman con¬ 
tinuum radiation. 

Based on these results, we rank each of our galaxies 
in each of the three above parameters from the galaxy 
having the strongest indication of escaping ionizing radi¬ 
ation to the one with weakest indicator. Thus, galaxies 
with a large residual flux, positive Reqw, and a large dis¬ 
tance from the star forming ridge in log ([S II]/Ha) vs. 
log([Oiii]/iL/?) space will have the first rank in each cat¬ 
egory. Then we average the rank in all three categories 
to order our LBAs from most (1) to least (22) likely to 
show evidence for escaping ionizing radiation. In one of 
our LBAs (J0921) - ranked “third” in our scale of indirect 
Lyman continuum escape indicators - follow-up observa¬ 
tions with HST COS have indeed revealed a dust-free 
escape fract ion of 21%, the highest measured in any lo¬ 
cal galaxy (Borthakur et al. 2014). Thus, our indirect 
indicators do (at least in this case) correctly predict the 
presence of real escaping ionizing flux. The differences in 
key absorption features for our “leaky” and “non-leaky” 
candidates are best shown in figure where we have 
combined the spectra of the “leakiest’Tialf of the sample 
as well as the least “leaky” half of the sample. It is easy 
to see evidence of the parameters we selected on, e.g. the 
“leaky” sample shows much weaker Si H absorption at 
1260 A but we also see evidence of other emergent prop¬ 
erties e.g. the equivalent width of Lyman a is smaller in 
the less “leaky” galaxies. 

5.2. Comparison with Galaxy Properties 

Given these results, we now compare our “leaky indi¬ 
cator” rank to other galaxy properties to look for those 
that correlate with the likely escape of ionizing radiation. 
For a comparison of the correlation between individual 
indirect indicators of escaping Lyman continuum radia¬ 
tion and galaxy properties see Table 

We begin by evaluating three other parameters that 
have been proposed to be indirect indicators of leaki¬ 
ness. 009 and Hll suggested that the ratio of the SFR 
derived from dust-corrected FUV luminosity to that de¬ 
rived from the dust-corrected Ha luminosity was such an 
indicator. Jaskot & Oey (2013) proposed using the flux 
ratio of the [OHI]A 5007/[OH]A3727 emi s sion lines, and 


nume rous authors (Jaskot & Oey 2014 Shibuya et al. 


2014) have proposed using the equivalent width of the 


Table 1 

Pearson correlation between escape indicators and 
general LBA characteristics (significance of deviation 
from zero). 


Indicator 

EWi 

SFRjrijv 2 
area 


T? 4 

c res 

0.42(0.026) 

0.28(0.153) 

0.41(0.031) 

R ^ 

J-^eqw 

0.42(0.021) 

0.22(0.142) 

0.20(0.171) 


0.18(0.236) 

0.77(0.001) 

0.63(0.001) 


^ Lyman a equivalent width (EW) measured over the same 
region as Reqw but calculated using only emission. More 
details can be found in section l3.5l 

^ Star formation rate (Mq yr“^) calculated from the FUV lu¬ 
minosity cor rected with the observed f ar-infrared luminosity 
according t o|Kennicutt Sz Evans] j2012[ For more information 
see section |4.1[ l divided by the area of the galaxy inside the 
petrosian radius enclosing 50% of object flux in the UV (kpc; 
for additional details see section |3.6[ l 

^ Outflow velocity calculated as the fir st ni oment of the 
Si III absorption line (km s“^; see section |3.2| for additional 
information). 

Residual Flux in the Si ii absorption line measured at the 
centroid of the Si ii 1260 A line (%). For a description of the 
measurement see section l3.3l 

^ Lyman a equivalent width ratio of blueshifted/redshifted 
equivalent width where a negative value is absorption and 
a positive value is emission. More details are provided in 
section fsTsI 

® Perpendicular distance from the star-forming ridge line 
measured on a plot of log (S ii/Ha) vs. log([OIII]//f,d)(dex). 

For additional details refer to section lT3l 

Figure We see no significant correlation of rank with 
either me SFR ratio or [OHI]/[OH]. We believe in this 
first case this is probably because the intrinsic scatter in 
these SFR indicators is higher than the escape fraction 
in these galaxies. On the other hand, there is a signif¬ 
icant correlation between rank and Lyman a equivalent 
width, albeit with exceptions like J0921 which has rela¬ 
tively weak Lyman a emission. 

Next, we examine how the leakiness rank correlates 
with the basic physical and dynamical properties of the 
galaxies. In order of the strength of the correlation, we 
find that galaxies which are most likely to be “leaky” 
tend to have significantly higher star formation rates per 
unit and higher measured outflow speeds (see Fig¬ 
ure [13). We do not find significant correlations between 


leakiness rank and either the SFR or the burst mass. 

6. DISCUSSION 

6 .1. Relation Between Lyman a Emission and Eseaping 
Ionizing Radiation 

There have been increased efforts to model the ex¬ 
pected Lyman a profile shape for different possible 
geometries and kinematic models for the interstellar 
medium. Much of the motivation has been to attempt 
to relate the properties of the Lyman a emission line to 
the esc ape of ionizing radia tion. Simpler models , such as 


Erb et al. (2010) already 
e origins of our Lyman 


Lyman a emission line as a diagnostic. We show the 
correlations of each of these with our “leakiness rank” in 


those of Steidel et al. (2010) and 
provide insight into the possib. 
a line profiles. A classical P-Cygni profile is the product 
of a spherical outflow suc h as those produced in the winds 
of ea rly-type stars (e.g. |Castor||1 970[ [Castor fc Lamers 
19791) or starburst galaxies (e.g. [Ahn et al.||2000| IVer- 
hamme et al.[[2006 ). We see blueshitted absorption from 
outflowing gas approaching the observer on the front side 
of the galaxy. Meanwhile, Lyman a photons are able to 
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Figure 10. Comparison of our indirect indicators of escaping Lyman continuum flux. Top, left: comparison of residual flux in 
the Si II absorption line (%) to the net equivalent width of Lyman a (ratio of equivalent width blueward:redward of the Lyman 
a rest wavelength). Top, right: comparison of residual flux in the Si II absorption line (%) to the distance from the star-forming 
ridge line in log (S Il/Ha) vs. log([OlIl]/iL(3) space (dex). Bottom, left: comparison of the net equivalent width of Lyman 
a (ratio of equivalent width blueward:redward of the Lyman a rest wavelength) to the distance from the star-f orming ridge line 
in log (S Il/Ha) vs. log([OlIl]/lT/3) space (dex). J0921, our object with confirmed Lyman continuum emission (|Borthakur et al.| 
20141, is represented as a yellow diamond. Objects for which a local continuum fit to Lyman a was used are represented here 
as purple squares. The Pearson r correlation coefficient is given in the top left corner of each figure with statistical significance 
of the deviation from zero in parentheses. 


resonantly scatter into our line of sight and escape the 
galaxy most effectively by scattering off outflowing ma¬ 
terial moving away from the observer on the back-side 
side of the outflow. These redshifted photons can then 
successfully traverse the HI on the front-side, which we 
see as redshifted emission. A low column density of neu¬ 
tral hydrogen gas in the galaxy (and thus a small r^ya 
at the systematic velocity of the galaxy and in outflows) 
in a spherically symmetric outflow would allow Lyman 
a photons to escape through an optically-thin medium 
even in outflowing gas on the side of the galaxy near the 
observer and produ ce blueshif t ed em ission instead. This 
idea was applied by Erb et al. (2010) to Q2343-BX418, a 
young. 


y galaxy at z = 2.3, which 
displays significant blueshifted emission in the Lyman 
a line. 

Much effort has been put into creating more compli¬ 
cated radiative transfer models of L yman a photons in an 


expanding m e dium . Most recently, Behrens et al. 


Duval et al. (12014); Verhamme et al. (2014) hn 


2014); 

that, 


tor an expanding shell with holes in the neutral gas, a et al.|20i4 ). Given the prorninence of the Lyman a. emis- 


Lyman a profile with a main peak at the systematic red- 
shift of the object as well as emission blueward of the 
systematic redshift are good indicators of escaping ioniz¬ 
ing continuum radiation- both characteristics seen in our 
Lyman a profiles where we predict a high escape fraction 
of ionizing radiation. While symmetric Lyman a pro¬ 
files coul d also be due to radiative transfer i n a static 
medium (Neufeld|1991 Verhamme et al.|2006 e.g.) this 
possibility is ruled out by the high outflow speeds we 
measure in the interstellar absorption lines. Most radia¬ 
tive transfer models show a secondary, redshifted peak 
in Lyman a from emission in an expanding shell which 
we see blended with the central peak of photons escaping 
directly into our line of sight from the clumpy ISM. 

Most importantly perhaps, we do find a statistically 
significant correlation between our three indirect indi¬ 
cators of escaping ionizing radiation and the equivalent 
width of the Lyman a emission line. A high equivalent 
width is often used as a lin e of evidence for escaping ion¬ 
izing conti nuum radiation (Jaskot & Oey||2014 Shibuya 
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Figure 11. Coadd of the 11 most “leaky” objects (black) and the 11 least “leaky” objects (red). No weighting of individual 
spectra was applied. All spectra have been moved to the same wavelength grid as the SB99 models and the best-fit stellar model 
subtracted. Important lines are noted. 


sion line and the many on-going and planned surveys 
of Lyman a emitters at high-redshift, this is result is a 
promising one for understanding re-ionization. 

6 .2. The Relation Between Starhurst Compactness and 
Escaping Ionizing Radiation 

Hll argued that LB As with a dominant, compact star- 
burst were the best candidates for high escape fractions 
of ionizing continuum radiation. They argued that these 
objects provide extreme feedback (in the form of radia¬ 
tion pressure, wind ram pressure, and a high intensity 
of ionizing radiation), and this allows them to blow out, 
ionize or otherwise disrupt the neutral ISM. We have 
confirmed and quantified that result, finding that our 
galaxies ranked as the strongest candidates for high es¬ 
cape fractions of ionizing continuum radiation have the 
highest rates of star formation per unit area though there 
is a less clear link between those objects officially clas¬ 
sified as DCOs and those that are “leaky”. While all of 
our galaxies are chosen to have compact star-forming re¬ 
gions, i?5o 1.0 kpc, it is our most compact sources that 
make the best candidates for leaking Lyman continuum. 
We also confirm the result in Hll that the leaky galax¬ 
ies are characterized by high wind outflow speeds. We 
show the connection between starburst compactness and 
outflow speed explicitly in Figure Thus, all of these 
relations create a consistent picture where concentrated 
massive starbursts are able to drive strong winds capa¬ 
ble of creating holes in the optically-thick ISM through 
which ionizing radiation might escape. 

7. CONCLUSION 

In this paper we have reported on our analysis of a 
sample of 22 local “Lyman Break Analogs” (LBAs) us¬ 
ing HST COS UV-spectroscopy plus ancillary ultraviolet, 
optical, and infrared data. Our goals were to use these 
data to better understand how to recognize galaxies that 
are likely to have escaping Lyman continuum radiation, 
and to obtain new insights into the conditions and pro¬ 
cesses that enable ionizing radiation to escape. We have 


argued that these goals can be best accomplished using 
relatively nearby galaxies. We have fit SB99 models of 
instantaneous and continuous starbursts to characterize 
the stellar population of each LBA, and find these are 
young starbursts with the prominent populations of O 
stars needed to produce substantial amounts of ioniz¬ 
ing radiation. We have subtracted the resulting model 
spectra of the starburst to create a spectrum that shows 
only features due to the interstellar medium (ISM). By 
examining the low-ionization ISM lines we see evidence 
for only partial coverage of the starburst by outflowing 
neutral gas: while the absorption features are saturated 
(optically thick), the flux in the line core does not go to 
zero. This is corroborated by our measurements show¬ 
ing cases in which a significant relative amount of the 
Lyman a emission line profile is blue-shifted relative to 
the galaxy systemic velocity. This indicates, most likely, 
that HI only partly covers the front-side of the galactic 
outflow. Finally, we use SDSS optical spectra to show 
that the [SH] emission line doublet is unusually weak in 
many of our galaxies (compared to typical SDSS galax¬ 
ies). This may indicate the presence of matter-bounded 
HII regions (gas that is optically thin to the Lyman con¬ 
tinuum) . 

We have found that these three indirect indicators 
of “leakiness” all correlate well with one another. We 
have therefore combined all three diagnostics to rank- 
order our galaxies in terms of likely leakiness. We have 
noted that the inference of leakiness is corroborated by 
the direct detection of escaping ionizing radiation from 
the galaxy we rank number 3 out of 22 in our sam¬ 
ple (Borthakur et al. 2014). We have also found that 
our leakiness rank is well-correlated with the equivalent 
width of the Lyman a emission line, which has been 
widely suggested as a diagnostic of escaping ionizing radi¬ 
ation. Two other proposed diagnostics did not correlate 
with our leakiness rank: the ratio of the star-formation 
rates measured using the far-UV vs. Ha dust-corrected 
(intrinsic) luminosities, and the flux ratio of the opti¬ 
cal [OHI]5007/[OH]3727 emission lines. Evidently other 
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Figure 12. Correlation between our leaky galaxy candidate ranking and various properties expected to correlate with leaky 
Lyman continuum. The top left of each figure lists the calculated Spearman rank correlation coefficient for each relation as 
well as the statistical deviation from zero (in parentheses). J0921, our object with confirmed Lyman continuum emission 
(Borthakur et al.|2014|, is represented as a yellow diamond. Objects for which a local continuum fit to Lyman a was used are 
represented here as purple squares. 1 (top, left): Ratio of SFR derived from IR-corrected FUV luminosity to extinction-corrected 
Ha luminosity. We find only a weak correlation between this value and our “leakiness” ranking. 2 (top, right): The equivalent 
width of Lyman a (measured in emission only) compared to our leaky galaxy rank. We find a strong correlation between the 
Lyman a equivalent width and our ranking. 4 (bottom, left): The ratio of [Olll] A5007,4959 flux to [Oil] A3727 emission to 
our leaky candidate ranking. A deficit of [Oll]is often used to identify potential leaking galaxies but we find an anti correlation 
between this indicator and our leaky candidate ranking. 


factors are more important in determining these ratios 
in our sample of galaxies. 

We then correlated our leakiness rank with the basic 
properties of the galaxy and its starburst. We found 
that the strongest correlations of leakiness were with the 
compactness of the starburst (SFR/area) and with the 
outflow speed measured using the ISM absorption lines. 
We also showed that the outflow speed is strongly corre¬ 
lated with the SFR/area. We have therefore argued that 
the extreme feedback associated with a compact star- 
burst is responsible for creating holes in the neutral ISM 
that allow ionizing Lyman continuum photons to escape. 
This is likely to be accomplished through a combination 
of the high intensity of the ionizing radiation field and 
the strong outward pressure exerted by the radiation and 
hot starburst-driven wind. 

A similar, but more ubiquitous population of compact 
star-forming galaxies at high redshift could create escape 
fractions high enough to re-ionize the universe. In addi¬ 
tion to providing these new clues to the processes that 


allow ionizing radiation to escape gas-rich galaxies, our 
analysis has established several indirect indicators that 
can be used to identify plausible “leaky” galaxies at high- 
redshift. This is crucial because the opacity of the inter- 
galactic medium makes it impossible to directly observe 
the escape of ionizing radiation from galaxies at or near 
the epoch of re-ionization. 
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